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Abstract 

The results of highly accurate calculations are presented for all twenty-two known bound S{L = 
0) — ,P{L = 1) — ,D(L = 2)— and F{L = 3)— states in the six three-body muonic molecular ions 
pp^,pdfi,ptfi, ddnjdtfi and ttfj,. A number of bound state properties of these muonic molecular 
ions have been determined numerically to high accuracy. The dependence of the total energies of 
these muonic molecules upon particle masses is considered. We also discuss the current status of 
muon-catalysis of nuclear fusion reactions. 

PACS number(s): 36.10.+Di, 36.10.-k and 31.10.-Fz. 



* E-mail address: lafrolov@uwo . cal 
'''E-mail address: ldwardlaw@uwo . cal 



I. INTRODUCTION 



In this study we consider the bound state spectra of the muonic molecular ions 
ppfi,pdij,,ptfi, ddn^dtfi and ttji. In this paper the notations p,d,t designate the nuclei of 
hydrogen isotopes (protium, deuterium and tritium, respectively), while /z means the neg- 
atively charged muon /i^. Our main goal is to determine total energies and other bound 
state properties in these muonic molecular ions to high enough numerical accuracy to be 
sufficient for all current and anticipated future experimental needs. 

In general, the bound state spectra in the six muonic molecular ions ppfi,pdfi,ptiJ,, 
ddjji, dtfi and tt/i can be separated into three different groups [2| on qualitative grounds. 
The first group includes the three light muonic molecular ions ppfi,pdij, and ptfi. Each of 
these systems has two bound states: one S{L = 0)— state and one P{L = 1)— state, where the 
notation L means the total angular momentum of the three-body system. Neither of these 
two states is weakly bound. Note that each of these light muonic molecular ions contains 
at least one protium nucleus. The second group includes the two 'intermediate' muonic 
molecular ions ddfi and dtfi each of which has five bound states: two S{L = 0)— states, 
two P{L = 1)— states and one D{L = 2)— state. One of these five states (the excited 
P*{L = 1)— state) in each of these ions is weakly bound. The third group contains only the 
heaviest muonic molecular ion ttfj, which has six bound states (and no weakly bound states): 
two S{L = 0)— states, two P{L = 1)— states, one D{L = 2)— state and one F{L = 3)— state. 

It can be shown (see, e.g., ^]) that the total number of bound states in any muonic 
molecular ion a^h^ ii~ is determined by the lightest nucleus in this ion. This explains why 
only three groups of different bound state spectra can be found in six such ions. Moreover, 
it follows that there must be a similarity between the energy spectra of the 'protium' muonic 
molecular ions ppii,pdfi and ptfi. An analogous similarity can be found in the bound state 
spectra of the ddfj, and dtfi ions. It can be shown that in such 'families' of muonic molecular 
ions the symmetric ion, e.g., ppfi, always has the maximal binding energy of the three 
protium ions: ppfi,pd^ and ptfi. By using these similarities between the bound state spectra 
in each of these 'families', one also finds a number of useful relations for the total and binding 
energies as well as for other bound state properties of different muonic molecular ions. For 
instance, let us assume that we know that the excited P— state in the ddfi ion is weakly 
bound and its binding energy is ~ -1.9745 eV (see, e.g., ^]). From the similarity of the 
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bound state spectra of the ddfj. and dtfi ions it one predicts that the corresponding excited 
P— state in the dtfi ion is also weakly bound and its binding energy is above -1.9745 eV, i.e. 
-1.9745 eV < e{dtfi, P*{L = 1)) < 0. This prediction turns out to be correct, as the dtfj, ion 
has a weakly bound excited P— state with the binding energy ^ -0.660 eV. 

Our labels for the bound states in muonic molecular ions are based on atomic LS- 
notations (see, e.g., 4!]). Note that there is another classification scheme which is still 
in use for muonic molecular ions and which was originally introduced to classify bound state 
spectra in adiabatic molecular ions, e.g., in the molecular ion j^. In this scheme each 
bound state is designated by its rotational J and vibrational u quantum numbers, i.e., we 
have the (J, z/)-states. The ground state in any muonic molecular ion is designated as the 
(0,0)-state, while the excited P— state in this scheme is denoted as the (l,l)-state, etc. Each 
of these classification schemes has its own advantages and disadvantages in applications to 
actual systems. Note also that there is a uniform correspondence between the 'atomic' and 
'molecular' classification schemes. 



II. THE HAMILTONIAN AND WAVE FUNCTIONS 

As mentioned above in this study we consider the bound state spectra in the six muonic 
molecular ions ppfi, pdfi, ptfi, ddfi, dtji and ttfi. All particles which form such three-body ions 
are assumed to be point and structureless. Each of these three particles has a finite mass 
which equals one of the masses m^,mp,md and/or rrit] the electric charges are = —1 and 
Qp = Qd = Qt = +^ (in muon- atomic units, where ^ = 1, = 1, e = 1). The Hamiltonian H 
of the three-body muonic molecular ions is written in the form, e.g., for the a^b^fi^ ion 

2m^Vm„ rrif, ^' Tab r^/, 

where Vj = (^-^r-, ^) ^) and i = a,b, /i. In muon-atomic units we have 

H = -U—^l + + + — (2) 

2Vma rrih Tab r^^, 

where the two masses and nif, of nuclei of two hydrogen isotopes must be expressed in 
terms of the muon mass m^. In fact, in this study only the muon-atomic units {h = 1, = 
1, e = 1) are used. Advantages of these units are discussed in Section IV below. 

Our computational goal is to determine exceptionally accurate solutions, i.e., the 
eigenstates and corresponding wave functions of the non-relativistic Schrodinger equation 



1, r2, vs) = E ■ \I/(ri, r2, r^), where E < and the non-relativistic Hamiltonian is writ- 
ten in the form of Eq.([2]). In actual calculations the wave functions of muonic molecular 
ions are usually approximated with the use of different variational expansions. In this work 
we shall consider the exponential variational expansion in relative coordinates ri2,ri3,r23 
6|. Here and everywhere below in this work the notation =\ r j — |= r^j designates the 
relative coordinate between particles i and j. In many cases, however, it is very convenient 
to introduce three new variables ui,U2,us which are called perimetric coordinates. They 
are simply related to the three relative coordinates: Ui = |(rjfc + — rjk), and therefore, 
Tij = Ui + Uj, where {i,j,k) = (1,2,3). The perimetric coordinates are truly independent, 
and each of them varies from to +oo. This significantly simplifies derivation of the explicit 
formulas for all matrix elements needed in highly accurate computations of the bound states. 
The explicit form of the exponential variational expansion in perimetric/relative coordinates 
is 

1 - ^ 

^LM = -(1 + /t^2l) C'i3^LAf (^31, r32)0i(r32, ^31, r2i) exp{-aiUi - (5iU2 - Ti^s) X (3) 

^ i=l h 

exp(2(5iMi + ieiU2 + ifiU^^) , 

where Ci are the linear (or variational) parameters, ai, /3i,ji,6i,ei and fi are the non-linear 
parameters and i is the imaginary unit. 

The functions 3^LM^(r3i, r32) in Eq-dS]) are the bipolar harmonics [3] of the two vectors 
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T31 = rsi ■ n3i and r32 = r32 ■ n32. The bipolar harmonics are defined as follows 

yi'^^\^,y) = x'^y'^ E Ci;l;..n^.>^^i™,(n.)y,,™,(n,) (4) 



where C|[mi-^2m2 Clebsch-Gordan coefficients (see, e.g., [3]) and the vectors ria, = ^ 

and Uy = ^ are the corresponding unit vectors constructed for arbitrary non-zero vectors x 
and y. Also, in this equation L is the total angular momentum of the three-body system, 
i.e. LF''^ LM = L{L + l)\E'iAf, while M is the eigenvalue of the operator, i.e. L^'^ lm = 
M^> LM- In actual calculations it is possible to use only those bipolar harmonics for which 
(■I + ('■2 = L + where e = or 1. The first choice of e (i.e. e = 0) corresponds to the 
natural spatial parity xp = of the wave functions. The second choice (i.e. e = 1) 

represents states with the unnatural spatial parity xp = {~'^)^^^- In this work we shall 
consider only the bound states of natural parity, since only such states exist in real physical 
systems. An additional family of polynomial-type functions ^1(^32, rsi, r2i) are also used in 



Eq.(|3]) to represent the inter-particle correlations at short distances. In general, these simple 
polynomial functions allow one to increase the overall flexibility of the variational expansion 
Eq.(|n]). In our present calculations, however, these additional functions were chosen in the 
form 0i(r32, Tsi, r2i) = 1 for i = 1, . . . , N. The operator P21 in Eq-© is the permutation of 
the identical particles in symmetric three-body systems, where n = ±1, otherwise k = 0. 
In this study we assume that k = xp = (^1)^ symmetric muonic molecular ions 

ppfi, ddfi, ttfi and k = for all non-symmetric ions pdfi,ptij,, dtfi. 

In general, highly accurate computations of bound states in muonic molecular ions are 
not easy to perform, since there are bound states with different angular momenta L [L 
= 0, 1, 2 and 3) and some of these states are very weakly bound. Variational expansions 
used in highly accurate calculations must provide fast convergence rates for each of the 
bound states, including all weakly bound states. Note that the actual goal of many current 
calculations of muonic molecular ions is the computation of various bound state properties, 
rather than the energies. In general, the convergence rate for some of these properties, 
including many nuclear- nuclear expectation values, e.g., the expectation values which include 
the nuclear-nuclear delta-function 5++ (e.g., the (5++^^—) expectation values for n > 1), is 
substantially slower than for the total energies. This explains our current need for highly 
accurate wave functions. For instance, as follows from numerical calculations to determine 
the ) expectation values, to the accuracy ±1 ■ 10^® one needs to use wave functions 

which provide an accuracy ~ 1 • 10^^^ a.u. for the total energy. Such values are needed 
in computations of the lowest order relativistic and QED corrections. A separate, but 
serious problem for accurate computations of three-body systems is the adiabatic divergence 
described in Isj. This problem always appears when 'pure atomic' variational expansions 
are applied to the two-center Coulomb systems and/or to the systems close to them. For 
our systems this means that the convergence rates observed for the ppfi and pdfi ions are 
relatively high in comparison to the analogous convergence rates for the dtfi and ttfi ions 
which are substantially slower. 

Nevertheless, variational calculations of bound states in muonic molecular ions are of 
great interest for predicting the physics of few-body systems as well as in some applications. 
In general, the study of bound state spectra in muonic molecular ions has provided us with 
a large amount of very valuable information and drastically improved our knowledge about 
the bound state spectra in arbitrary Coulomb three-body systems. Furthermore, all muonic 
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molecular ions are three-body systems with unit charges. The energy spectra in such systems 
have many significant differences from known atomic spectra. In particular, any Coulomb 
three-body system with unit charges has a finite number of bound states The only 
exception to this rule is the °°H^ ion which has infinite number of bound states j^. 



III. OPTIMIZATION OF THE NON-LINEAR PARAMETERS IN TRIAL WAVE 
FUNCTIONS 

Optimization of the non-linear parameters in trial wave functions is a central part of 
construction of any highly accurate solution of the Schrodinger equation. By optimizing the 
non-linear parameters in some rapidly convergent variational expansions one can produce 
extremely accurate and compact wave functions for arbitrary three-body systems. Such 
wave functions can be used in further calculations of various bound state properties and 
different corrections to the non-relativistic energies of three-body systems. In this study we 
construct our highly accurate variational wave fuuct.ons. E,.©, wth the use of a two-stage 
optimization strategy of the non-linear parameters a]. The two-stage procedure consists of: 
(a) construction of a short-term booster wave function with carefully optimized non-linear 
parameters, and (b) quasi-random choice of the remaining non-linear parameters from a few 
(usually, three - five) optimal boxes (or parallelotops) . To simplify the description of our 
optimization procedure below we shall assume that all exponents in Eq.(l3]) are real. In other 
words, all parameters 6i, Cj and fi in Eq.([3]) equal zero identically for i = 1, . . . , A^. In fact, 
these later parameters are really needed only for highly accurate calculations of the pure 
adiabatic three-body systems and systems close to them. Among all muonic molecular ions 
ppfi,pdij,,ptfi, ddij,,dtfi and ttfi even the two heaviest systems dt/i and ttfi are not close to 
adiabatic systems, examples of which are the one-electron DT+ and T^ ions. 

For muonic molecular ions considered in this study the short-term booster wave function, 
Eq-dS]), usually includes Nq = 400 - 600 basis functions (exponents) with 3A^o = 1200 - 1800 
non-linear parameters in them. All these parameters must carefully be optimized. After 
such an optimization the short-term booster function provides 11-15 exact decimal digits 
for each bound state energy in the considered muonic molecular ions. It appears that the 
overall accuracy of such short-term wave functions is much better for the protium muonic 
molecular ions ppfi,pd^,ptn than for the heavier ions ttfi and dt^. At the second stage 
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of our optimization procedure the remaining 3{N — Nq) non-linear parameters in the wave 
function Eq.(|3]) are chosen quasi-randomly from three different boxes or parallelotops. The 
total energies and other bound state properties obtained with such trial wave functions 
depend upon the boundaries of these boxes. In reality the boundaries of these three boxes 
can be described 16] with the use of 28 non-linear parameters only. The numerical values 
of these 28 parameters wp:e optimized approximately with the use of = 800, 1000, 1200 
and 1400 basis functions [6] . These values allowed us to determine the approximate limit for 
each of these 28 parameters as — > oo by extrapolation. These limiting values have been 
used in our final computations. ^ 

The two-stage strategy proposed in [2] and described above allows one to obtain very ac- 
curate variational wave functions based on the use of exponents in relative and/or perimetric 
coordinates. In particular, the overall accuracy of our results obtained in this study for all 
considered muonic molecular ions (see Tables I, II, III and IV below) is significantly higher 
than the accuracy obtained for these ions in earlier calculations. Moreover, by using this 
optimization strategy we expect to be able to increase the accuracy in future calculations by 
a factor of ~ 10^ — 10^ which would be sufficient for all future anticipated theoretical needs. 
The described two-stage optimization procedure has been used in this study for all bound 
S{L = 0)— and P{L = 1)— states in muonic molecular ions, including all excited states. For 
the bound D{L = 2)— states in heavy muonic molecular ions {ddfi,dtfi and ttfi) we used 
another approach in which the short-term booster function is not constructed. However, 
it is clear that our current strategy for optimization of the non-linear parameters in the 
wave functions of the bound D{L = 2)— states is not optimal. In future studies we want to 
improve this strategy and produce highly accurate results for all bound D{L = 2)— states in 
muonic molecular ions. 



IV. VARIATIONAL ENERGIES OF MUONIC MOLECULAR IONS 

The results of our calculations of different bound states can be found in Tables I - VI. 
All computations are performed in muon-atomic units, where = l,h = 1 and e = 1 and 
we used the following values of nuclear masses [l^ , jilt : 

rrif, = 206.768262me , nip = 1836.152701me (5) 
rud = 3670.483014me , rut = 5496.92158me 
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where rrie designates the electron mass. Note that our highly accurate computations in this 



study are performed with the use of 84 - 104 decimal digits per computer word [12|, |l3| . 
allowing total energies to be determined to an accuracy ~ 1 ■ 10^^° — 1 ■ 10^^^ m.a.u. A 
natural and effective way to perform such calculations is to assume that all particle masses 
and corresponding conversion factors (e.g., the factor Ry below) are exact. Such assumptions 
are always made in papers on highly accurate computations in few-body systems (see, e.g., 
[l^ and 15|). The known experimental uncertainties in particle masses and conversion 
factors are taken into account at the last step of calculations, when the most accurate 
computations are simply repeated for a few times with the use of different particle masses 
and conversion factors. Analogously, the lowest order relativistic and QED corrections can be 
determined as the expectation values of some operators computed with our non-relativistic 
wave functions. To avoid a substantial loss of numerical accuracy during such computations 
these non-relativistic wave functions must be extremely accurate. 

Table I contains the total variational energies obtained for the ground S{L = 0)— states 
of the non-symmetric muonic molecular ions pdfi,ptfi and dtfi. Table II includes the total 
energies for the excited S*{L = 0)— states of the dd/j,, tt^ and dt/j, muonic molecular ions. The 
total energies of the ground S{L = 0)— and 'rotationally' excited P{L = 1)— states of the 
symmetric muonic molecular ions ppfi, ddfi and ttfi are presented in Table III - IV. Table V 
contains the best variational energies obtained in our computations of the D{L = 2)— states 
of the ddfi, ttfi and dtfi muonic molecular ions. Highly accurate computations of these bound 
D{L = 2)— states have been performed with the use of extended arithmetic (for discussion 
and references, see, la]). In earlier works the bound D{L = 2)— states have been determined 
with the use of quadruple precision only. However, our current total energies determined 
for the bound D{L = 2)— states are not as accurate as the total energies obtained for the 
bound S{L = 0)— and P{L = 1)— states (see comment above). 

Note also that the total energies of the P{L = 1)— states of the non-symmetric muonic 
molecular ions pd^.pt^ and dt^ and excited P*(L = 1)— states of the ddfi,ttfi and dtfi ions 
have recently been determined in 16|]. The most accurate variational energies for all known 
22 bound states in the set of six muonic molecular ions studied here (expressed in muon- 
atomic units) can be found in Table VI. Note that the F(L = 3)— state of the tta ion has 
not been re-calculated in this study, but instead it has been taken from our earlier work [2| , 
where this state was computed with the use of quadruple precision. 
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As follows from Tables I - VI the total energies obtained in this study for different bound 
states in six muonic molecular ions are significantly more accurate than the corresponding 
energies computed in earlier studies (see, e.g., j2l, [l4| and 15|). The current wave functions 
are more compact and have better overall quality than wave functions obtained in 



and 



15| . They can be used for highly accurate computations of other bound state properties. 



including properties which contain singular expectation values. 

Our variational wave functions can be used to compute some bound state properties of 
muonic molecular ions. A large number of bound state properties have been computed 
in our earlier studies (see, e.g., [sj] and references therein). However, some of the bound 
state properties could not be determined to high numerical accuracy, due to relatively low 
accuracy of the wave functions used in earlier studies. It was clear that the expectation 
values of some nuclear- nuclear properties, e.g., all properties which include the nuclear- 
nuclear delta-function, needed to be re-calculated with more accurate wave functions. By 
using highly accurate wave functions obtained in this work we have performed numerical 
re-calculation of a number of bound state properties for different muonic molecular ions. 

The computed expectation values can be found in Tables VII and VIII. Results presented 
in Table VII illustrate convergence of some expectation values upon the total number N of 
basis functions used. By comparing the expectation values computed with the use of different 
number(s) of basis functions we have determined the corresponding asymptotic values for 

— 7- oo. Formally, this procedure allows one to determine the number of stable decimal 
digits for each of the computed expectation values. Table VIII contains some bound state 
properties determined for the ground S{L = 0)— states in the pdfi and dt/j, muonic molecular 
ions and for the excited S{L = 0)— state in the dtfi ion. 



V. MASS DEPENDENCE OF THE TOTAL ENERGIES FOR MUONIC MOLEC- 
ULAR IONS 

Let us consider the problem of mass dependence for the total energies of muonic molecular 
ions. Briefly, this problem is formulated as follows. The Hamiltonian H of any muonic 
molecular ion, Eq.(Il]), contains three different masses ma^mbim^. These particle masses 
are the subject of constant experimental revision. Furthermore, various authors often use 
slightly different particle masses in their calculations. For highly accurate computations this 
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means the almost constant necessity of re-calculation of the energies and corresponding wave 
functions. In fact, all masses of particles which form muonic molecular ions are currently 
known to relatively high experimental accuracy, i.e. all possible 'mass-corrections' must be 
very small. This means that such corrections can be considered by using various methods of 
perturbation theory. Nevertheless, it is very interesting in some cases to evaluate corrections 
produced by the corresponding 'mass shift'. The most important bound state property is 
the total energy E. The first order variation of E with the particle masses is written in the 
form 

E = E -\- fy 

-'-'new -'-'our ~ 



■ m^jnew) m^(oMr) j ^ ^r m^(new) m^(oMr) j 
^nitinew) mt{our) ^ '-rridinew) md^our)^ 

in the case of the dtfi ion. In Eq.(l6]) and Eq.(I7]) below all energies must be expressed in 
the same units, e.g., in atomic units or in muon-atomic units reduced to the same muon 
mass. The same formula can be written for any non-symmetric muonic molecular ion. In 
such cases we always have a > /3, if the coefficient a corresponds to the mass shift produced 
by the heaviest nucleus. For symmetric muonic molecular ions, e.g., for the ddfi ion, the 
analogous formula is 



E = E -I- o- 

-'-'new -'-'our ~ " 



(7) 

-md[new) ma[our)-' 

In these equations the notation 'our' denotes the mass value used in this study, while the 
notation 'new' designates a different mass value, e.g. from some work performed in the 
future. In general, the numerical values of the coefficients a and /3 in Eq. ([6]) (also called the 
mass gradients) are determined from separate energy calculations with different masses. In 
our earlier works we have used a very simple approach based on four additional calculations 
with different masses for non-symmetric muonic molecular ions [l^ . For symmetric muonic 
molecular ions one has to perform at least two additional calculations with the two different 
mass ratios. This method is simple, but it is not very accurate. 

Recently, we have developed a more accurate procedure. To describe this procedure 
let us consider the ground P(L = 1)— state in the ddn ion. The mass ratio "^1^^°'^'^}^ is 

<=> \ / r~ md(our) 

designated below as xn, while the notation h stands for the difference _ 

° m^{new) mc^(our) 

used in Eq.(I7j). In these notations the total energies Eour and Enew are Eour = E{xo) and 
Enew = E{xo -|- h), respectively. The absolute value of h is assumed to be small. Our new 
method for calculation of the mass gradient a is based on calculation of the four ground 
state energies E{xq — 2h), E{xq — h), E{xq -^h) = Enew and E{xq -|- 2h). These four energies 
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are determined in calculations performed with the use of the maximal number of basis 
functions (at the maximal dimension). No re-optimization of non-linear parameters in the 
wave function is required during such additional computations. The 'mass gradient' a in 
Eq.(I7]) can now be computed with the use of the formula 



a 



12h 



1 

E{xo-2h)-8E{xo-h) + 8E{xo + h)-E{xo + 2h) + —A (8) 



where v4 is a numerical parameter of order of unity. It can be shown that this parameter 
equals the fifth order derivative of the total energy E with respect to the mass ratio ^ 
computed in some point between xq — 2h and Xq + 2h. The formula, Eq.®, has a very good 
numerical accuracy in actual applications to all bound states, except weakly-bound states 
(discussed below). 

The results of numerical calculations of the mass gradients a in Eq.(IH]) for the P{L = 
1)— states in the symmetric ppfi,ddfi and ttfi muonic molecular ions are presented in Table 
IX. Table IX also contains the results for the excited P*(L = 1)— state in the ttfi ion. All 
these energies have been re-calculated to the 'standard' muon mass = 206.768262 rUe- 
This Table also contains all intermediate energies E{xq — 2h),E{xQ — h),E{xo + h) and 
E{xq + h) needed in such calculations. As follows from Eq.® our mass gradients a are 
accurate to 11 - 13 decimal digits. The same procedure can be used for other bound states 
in muonic molecular ions. 



A. Mass dependence for the weakly-bound states 

The mass dependence of the total and binding energies for weakly bound states cannot 
be investigated with the use of the method described above, since the total energies E{xo + 
h) and E{xq + 2h) needed in Eq.([8]) can correspond to unbound states. This indicates a 
necessity to study the mass dependence of the total and binding energies for weakly bound 
states. Formally, to answer this question one needs to perform a number of calculations with 
different particle masses. For simplicity, let us discuss the weakly bound excited P*{L = 
1)— state in the dtfi ion. If the muon mass and the muon-deuterium mass ratio t = ^ 
increase, then the P*{L = 1)— state in the dt/j, ion will be less and less bound. Finally, at 
some critical values of and r the P*{L = 1)— state in the dtfi ion will become unbound. 

The approach based on the use of 'conventional wisdom' allows one quickly to find the 
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numerical value of the critical muonic mass. Indeed, the Hamiltonian of muonic molecular 
ions, Eq.(|2]), is a linear function of the muonic mass and/or the mass ratio t = 
Therefore, by performing calculations of the total energy of some bound state for different 
muonic masses one finds the dependence E{m^), or E{mr). The next step is to determine 
the critical mass for which E{m^) = , i-e. the binding energy of this state equals 

zero. For a very short interval of variation of the function E{m^) is almost a linear 
function upon the muonic mass m^, then it is relatively easy to obtain an approximate value 
of the critical muonic mass. 

This approach was used in some earlier works, but it produces incorrect results in ap- 
plications to the bound states which can disappear (as bound states) during variations of 
some physical parameter, e.g., the mass of the particle. For such states one needs to use an 
alternative approach described in [l8|. In this method the explicit expression for the total 
energy of the two-body 'almost unbound' system takes the form 

32 m^m? ^ ' 

where U = is the effective two-body potential in the two-body dji + d system. The 
value Ui corresponds to the mass mi at which e: = 0. The difference of masses nii and 
m is assumed to be small in comparison with each of the two masses (mi and m). The 
formula, Eq.(|9]), predicts a quadratic dependence of the energy e upon the mass difference 
A = nil — m, i-e. e{mi,m) = e{mi,A) ^ a ■ Although the actual e:(mi, A) dependence 
must be described by the formula Eq.([n]) to good accuracy, in reality however it is difficult 
to obtain the quadratic dependence for the e{mi, A) function (on the parameter A) by using 
the results of numerical calculations with the same number of basis functions. It is clear that 
in accurate calculations of bound states with A — one needs to increase the total number 
of basis functions and constantly re-optimize the non-linear parameters of the method. If 
such conditions are obeyed, then we can observe the expected quadratic dependence of the 
binding energy e: on A. Otherwise, from the results of such calculations one can see only an 
approximate linear dependence. 

Note that for real muonic molecular ions the total energy of the two-body system e given 
in Eq.Q is, in fact, the binding energy of the three-body ion, e.g., the rft/i ion, which 
corresponds to the lowest-by-energy decay channel: {dtn)^ = t/x(ground state) + d'^. For 
the dta ion the mass m in Eq.(l9]) is the 'reduced' mass m = while mi = ("■t+^-A'Hrf 

r~ \<—v mt+m^+nid ' ^ mt+inij,+md 
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is the 'threshold' mass, i.e. the reduced muonic mass for which the binding energy equals 
zero. 

The results of our calculations for the P*(L = 1)— state in the ddjj, ion can be found in 
Table X. In these calculations we have increased (at each step) the mass of the /i^ muon 
(see above) by one electron mass mg. As follows from Table X the binding energy of the 
P*{L — 1)— state in the dd/i ion decreases to zero. The value of muon mass at which the 
corresponding bound state becomes unbounded is called the threshold mass. For the ddfx 
ion such a muon threshold mass is designated as 771^. In calculations for Table X we have 
used N — 3300 basis functions. The £(mi, A) dependence which follows from the results 
of these calculations can be represented as a linear function. The threshold value of 
found from this linear dependence approximately equals to the value rh^ = + 7.95me ~ 
214.718262me. As follows from Table X the binding energy of the P*{L — 1)— state in the 
ddp, ion is -0.00002551843452 eV. This energy was obtained with N — 3300 basis functions. 
However, if we take 3840 basis functions, then the binding energy of the P*(L — 1)— state 
in the same ddfj, ion is fa -0.00007865184962 eV. i.e., it corresponds to the state which 
is slightly better bound. The muon threshold mass is always shifted to larger values of 
771^ if the total number of basis functions increases. This example explains the occurence 
of what is often called "running" mass threshold, a non-physical effect known to exist in 
numerical calculations of any weakly-bound state. It is interesting to note that the muon 
mass — 214.718262me exceeds our 'standard' muon mass (206.768262 nie) already by fa 
4 %, but the P*{L — 1)— state in the ddfj, ion is still bound (earlier estimations for such a 
mass deviation were around 0.5 - 1 %). 

Analogous results for the weakly bound P*{L = 1)— state in the rft/i ion can be found 
in Table XI. In computations performed for this Table we varied only the muon mass m^, 
while the deuterium and tritium masses have not been changed. At such conditions the muon 
threshold mass was found to be equal ~ m^^ -|- 1.99me. It is clear that the binding energy 
of the P*{L — 1)— state in the dt/j, ion is the function of the two independent mass ratios ^ 
and Therefore, Table XI gives only an approximate picture of how the total and binding 
energies (in eV) of the P*{L = 1)— state in the dtfi ion vary when the muon mass changes. To 
study the pre-threshold mass dependence of the weakly bound P*{L = 1)— state in the dtfi 
ion in detail one also needs to consider changes of the deuterium and tritium masses. Such 
calculations, however, are very difficult to perform, since they require substantial computer 
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resourses. 



VI. MUON STICKING PROBABILITIES 

Originally, all numerical computations of bound states in muonic molecular ions were 
motivated by various problems of muon-catalyzed nuclear fusion. In fact, the bound state 
computation of muonic molecular ions is only one of many problems which must be solved 
before we can discuss a possibility to use muon-catalyzed nuclear fusion for energy production 
and for other purposes. It is clear that the most interesting and promising case is the 
muon catalysis of the ((i, t)— nuclear reaction in the dtfi ion. A central problem here is to 
determine the muon sticking probability during the nuclear reaction dtfi = ^B.e+fi + n since 
the numerical value of this coefficient essentially determines the feasibility of using muon 
catalysis of nuclear fusion reactions for energy production purposes. Let us evaluate the 
muon sticking probabilities for this ion by assuming that the nuclear dt-fusion occurs only 
in the two S{L = 0)— states (ground and excited), ignoring the possibility of nuclear fusion 
in the P{L = 1), P*{L = 1) and D{L = 2) bound states of the dtfi ion. 

The analytical expression for the muon sticking probability for the bound S{L = 0)— state 



also 191 and 



20|) 



(initial state is designated with the subscript in; final state by fi) takes the form {g] (see 



r+oo 

Pin-ji = 47r(2£ + 1) I / <PU(^-'r)UQr)Rr,eirydr \^ (10) 
Jo 

where n and i are the appropriate principal and angular quantum numbers for the final 
hydrogen-like {'^B.e^)'^ ion with radial function Rne{r). The choice of the factor a in Eq. (fTO|) 



is discussed below. The ji{Qr) function is the spherical Bessel function (see e.g., |2l|): 



= ^^Ji^{x) (11) 



The factor Q is 



Q = m^v = J- — = 5.825011748 12) 

^ V (l + ^4)(l + "2„ + M4) ^ ^ 

where = 1838.683662 me is the neutron's mass, AE is the total energy release during 

the nuclear (d, t)— reaction and M4 = 7294.2296me is the mass of the ^He nucleus. 

In the formulas presented above (pmia'^r) is the initial 'post-process' wave function, i.e. 

the wave function of the system which arises when the sudden process (i.e. nuclear fusion) 
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is over. The function (f)in{a~^r) can be found from the bound state wave function of the 
initial three-body system. For instance, in the case of nuclear fusion in the S{L — 0)— state 
of the dt/i muonic molecular ion one finds: 



N 



4>in{rz2) = ^21 • *(r32, rai, rsi) = *(r32, rga, 0) = ^ Cj exp(-(Q!j + A)'^32) , (13) 

1=1 

where 821 is the nuclear delta-function and Cj are the linear variational parameters from 
Eq.(l). These coefficients have been determined during numerical solution of the Schrodinger 
equation for the initial three-body system (see Section IV). Note also, that after the 'sudden' 
nuclear fusion the new ^He nucleus arises at the same point '2'. This does not change 
the relative r^^ coordinate which is mass independent. After the nuclear reaction the 
coordinate becomes the helium-muonic relative coordinate. In two-body atomic problems, 
however, it is more convenient to use the mass-weighted coordinate r. The relation between 
the relative coordinate and mass- weighted r coordinate (which corresponds to the helium- 
muonic ion) is written in the form: 

m^M4 M4 
r = ar32 = -JT^^'^ , or r ^ r32 14 

where in muon-atomic units = 1 and M4 is the nuclear mass (in muon-atomic units) of 
the ^He nucleus. Finally, the initial wave function (j)in{a~^r) takes the form: 

N 

(l>in{(^~^r) =J2Ciew{-{<^i + I3i)a~^r) , (15) 
1=1 

where = For the muon and nuclear masses indicated above, one finds that 

a-\'^He) = 1.028346555. 

The formulas given above allow us to determine the muon sticking probabilities for the 
ground and excited S{L — 0)-states in the dt/j, ion. In these calculations we have used our 
best variational wave functions, as obtained in this study. After a number of calculations we 
have evaluated the total muon sticking probabilities Pg for the ground and excited S{L — 0)- 
states in the dtn ion. These two numerical values are close to each other. In fact, we have 
found that each of these two vahies (P^) is bounded in the interval 0.008923(3) < < 
0.008938(3) (for both bound S*— states). This means that the total number of rft— nuclear 
reactions catalyzed by one muon is k. — w 112. This means that before its decay 
one fj,~ muon can catalyze up to 112 nuclear reactions of c/i— fusion. This number is used 
in the next Section to discuss the possibility of exploiting muon catalyzed fusion for energy 
production purposes. 
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VII. CURRENT STATUS OF MUON-CATALYZED NUCLEAR FUSION 



Our evaluation of the factor k = ^112 performed above is based on an assumption 
that all nuclear (rf, t)— fusion reactions proceed from the bound S"— states of the dtiJ, muonic 
molecular ion. This factor does not take into account the stripping of muons from rapidly 
moving He nuclei and a possibility of nuclear (d, f)— fusion in the bound P— and D— states 
of the dtfi. Formally, an average time required to complete one (rf, t)— fusion reaction in 
any of the P— and/or D— states is substantially longer than the analogous time for the 
5"— state(s). In fact, such a 'reaction time' is comparable with the transition time which is 
needed to complete all transitions from the bound P{L = 1), P*{L = 1) and D{L = 2) states 
into the corresponding S{L = 0)— state(s). This directly follows from the fact that actual 
deuterium-tritium distances are larger in the P— or D— states of the dtfi ion than in the 
iS— states. Briefly, this means that all dtfi ions originally formed in the P— and D— states 
will emit radiation and make transitions into the corresponding S{L = 0)— states. In other 
words, a possibility of nuclear fusion from the 'rotationally excited' bound P— and D— states 
can be ignored and all reactions of the nuclear (rf, t)— fusion in the dt/j, ion can reasonably 
to assume to proceed only from the two bound 5*— states of the dtfi ion. Let P be exact 
sticking probability of the muon in the dtfi ion. The analogous value Pg is the muon sticking 
probability determined for the same dtfi ion, but only for its bound S{L = 0)— states. As 
follows from the discussion above we can replace the factor k = P"^ by the approximate 
value K ^ P^^. 

There are also a few other corrections which can change the numerical value of the factor 
K = P~^ ^ 112. The largest of such corrections corresponds to the 'muon stripping' during 
collisions of the fast (^He^)^ ion with neutral hydrogen molecules. However, even such a 
correction cannot change the predicted value of k by 40 % [1]. In other words, the maximal 
value of K is ~ 160 - 170. On the other hand, to reach break-even, i.e. to compensate for the 
energy spent for creation of one fi muon 8000 MeV j22|), one muon needs to catalyze 
at least 2285 nuclear dt— reactions. In this evaluation we have ignored all possible energy 
losses and assumed 100 % efficiency for each muon. In reality any thermal-to-electrical 
conversion has only ~ 30 % efficiency and only ~ 70 % of all muons can produce the 
maximal number of fusion reactions. With all these corrections one finds that the factor 
K must be ~ 8,000 - 11,000 to reach break-even. Such values are ~ 65 times larger than 
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the maximal value of k which has been measured experimentally (k ~ 150). If somehow in 
future experiments the numerical value of k will be increased up to 500, even then it will 
be ~ 20 times smaller the value which is needed to reach break-even. This indicates clearly 
that muon catalysis of nuclear reactions cannot be used for energy production purposes. 

It should be mentioned that originally the idea to use //"-muons for production of repet- 
itive nuclear reactions between light nuclei of hydrogen isotopes was proposed more than 
sixty years ago [23^. Based on an obvious chemical analogy these processes were called the 
muonic catalysis of nuclear reactions. It was confirmed in 2J] experimentally by observing 
two consecutive (p, ci)— nuclear reactions catalyzed by the same muon. The first numerical 
computations of the bound states in three-body muonic molecular ions were performed by 



Belyaev et al in 1959 



25| who found only 20 bound states in six ions ppfi,pdij,,ptfi, ddp, dtji 



25| was very low and the authors 



and tt/i. The overall accuracy of the procedure used in 
could not confirm the boundness of the excited P*{L = 1)— states (or (l,l)-states) in the 
ddfi and dtp ions. It was concluded only that, if such states are bound, then they are very 
weakly bound. The binding energy of these two states was expected to be smaller than 4.5 
eV, i.e. smaller than the binding energy of a typical molecule. Immediately after publication 



of 



25| an intense stream of speculations started about a possible interference (or resonance) 



between the formation of excited P*{L = 1)— states (or (l,l)-states) in the ddfi and dt/j, 
muonic molecular ions and different atomic/molecular processes in surrounding molecules 
(see, e.g., 26|] and references there in). Finally, in a few experimental studies performed by 
Bystritskii et al (see js^l and 28 1 and references therein) it was shown that one muon can 
catalyze approximately 10 - 20 {d, d)— nuclear reactions in liquid deuterium (D2) and 90 - 110 
(d, t)— reactions in the liquid equimolar deuterium-tritium mixture (D2 : T2 = 1:1). Such 
very large numbers of nuclear reactions catalyzed by one muon can be explained only by 
the resonance (or very fast) formation of dd/j, and dt/j, muonic molecular ions. Correspond- 
ingly, these processes were called 'resonance' muon- catalyzed fusion of nuclear reactions, in 
contrast with the 'regular' muon-catalyzed fusion observed in 24 1. 

In experiments performed in 1980's the total number of nuclear reactions catalyzed by one 
muon (i.e. the numerical value of the factor k defined above) for the equimolar deuterium- 
tritium mixture were evaluated a s ~ 150 (see discussion and references in [l|). This value 
is ~ 15 times smaller than the value which is needed for theoretical break-even and ^ 65 
times smaller than necessary for actual break-even. Therefore, we have to conclude that all 
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discussion of the 'bright future' for apphcations of the resonance muon-catahzed future for 
the energy production purposes appears to be groundless. 



VIII. CONCLUSION 

We have considered the problem of highly accurate calculations of bound states in the 
three-body muonic molecular ions ppfi,pdij,,ptfi, ddfj,, dtfi and ttfi. The study of bound state 
spectra in the muon-molecular ions is of interest for solving some theoretical problems and 
in a number of applications. In fact, our present knowledge of the bound state spectra 
in Coulomb three-body systems with unit charges is essentially based on knowledge of the 
spectra of the muonic molecular ions. Note that all muonic molecular ions can easily be 
created in real experiments and their various properties can be measured quite accurately. 
From a certain point of view, theoretical and experimental study of these ions is more 
interesting and informative than the traditional analysis of atomic three-body (i.e. two- 
electron) systems. 

The results of variational computations of the total energies for various bound states in 
the various muonic molecular ions are presented in Tables I - VI. Table VI contains the 
most accurately known predictions for variational energies in muon-atomic units for each 
of 22 known bound states for muonic molecular ions. The accuracy of these total bound 
state energies significantly exceeds the accuracy achieved in earlier studies. It should be 
mentioned that first variational calculations of muonic molecular ions started almost 45 



years ago 



29j, |30|, |3l|, |32j. In these works only S{L = 0)- and P{L = l)-states of 



muonic molecular ions were considered. It is interesting to note that at that time the non- 



variational calculations (see, e.g., 33| and references therein) of muonic molecular ions had 
a comparable overall accuracy for many bound S{L = 0)— and P{L = 1)— states in muonic 
molecular ions. In addition to this, the non-variational methods allow one to determine the 
total energies of the bound D{L = 2)— and F{L = 3)— states in muonic molecular ions 33|. 

Our first variational calculations of muonic molecular ions started 25 years ago 34 1. 
In particular, the first successful variational computations of the weakly bound P*(L = 



1)— states in the ddfj, and dtfi ions were performed in our work [34] and also in [35|. However, 
at that time we could not compute the bound D— and -F— states in the ddfi, dtfi and tt/j, ions. 
Furthermore, our maximal accuracy achieved at that time was relatively low. Currently, the 
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same energies for all S{L = 0)— and P{L = 1)— states in muonic molecular ions obtained 



m 



34| can be reproduced with the use of only 20 - 30 exponential basis functions in Eq.([3]), 
with carefully chosen non-linear parameters and ji in each basis function. The first 

variational computations of the bound D{L = 2)— states in the dt^^dd^ and tt^ ions were 
performed in 1986 fsG], while analogous calculations of the F{L = 3)— state in the ttfi ion 



were conducted 15 years later |2|. The bound D— state in the dtfj, ion was also calculated 
(variationally) by Kamimura in 1988 js?]. 

By using our highly accurate wave functions we have determined the expectation values 
of some bound state properties of muonic molecular ions. We also discuss the problem 
of mass shifts in muonic molecular ions, including the case of weakly bound states. A 
separate (but very important) problem is to evaluate the muon sticking probabilities for the 
S{L = 0)— states in the dtfi ion. By using these probabilities we estimated the total number 
of the reactions of {d, t)— fusion which can be catalized by one fi~ muon in a liquid equimolar 
D2:T2 mixture. The current status of the 'resonance' muon-catalyzed nuclear fusion is briefly 
discussed. It is shown that this process cannot be used for energy production purposes. 
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TABLE I: The total energies (E) of the ground bound S{L = 
muonic molecular ions in muon- atomic units (m^ = l,h = l,e = 
basis functions used in Eq.(2). 



0) — states in the non-symmetric 

1) . N designates the number of 



N E{pdfi) E{ptfi) E{dtfi) 

3300 -0.512 711 796 501 514 087 -0.519 880 089 782 919 741 -0.538 594 975 061 413 075 

3500 -0.512 711 796 501 514 175 -0.519 880 089 782 919 852 -0.538 594 975 061 413 084 

3700 -0.512 711 796 501 514 250 -0.519 880 089 782 919 605 -0.538 594 975 061 413 089 

3840 -0.512 711 796 501 514 283 -0.519 880 089 782 920 006 -0.538 594 975 061 413 095 

A -0.512 711 796 501 514 55(1) -0.519 880 089 782 921 0(5) -0.538 594 975 061 413 3(2) 

E"" -0.512 711 796 501 509 [2j -0.519 880 089 782 914 [2] -0.538 594 975 061 413 [2] 

E"" -0.512 711 796 500 8 [17] -0.519 880 089 781 9 \17] -0.538 594 750 58 flT] 

("^The best variational energies known from earlier calculations. 

TABLE XL The total energies (E) of the excited bound S{L = 0)— states in the muonic molecular 
ions in muon-atomic units (m^ = 1,^ = l,e = 1). designates the number of basis functions 
used in Eq.(2). 



N 


E{ddn) 


E{dtfi) 


E{ttn) 


3300 


-0.479 706 380 368 904 996 317 


-0.488 065 357 851 706 559 


-0.496 762 894 249 561 323 055 


3500 


-0.479 706 380 368 904 996 742 


-0.488 065 357 851 706 622 


-0.496 762 894 249 561 323 901 


3700 


-0.479 706 380 368 904 997 024 


-0.488 065 357 851 706 665 


-0.496 762 894 249 561 324 665 


3840 


-0.479 706 380 368 904 997 240 


-0.488 065 357 851 706 683 


-0.496 762 894 249 561 325 040 


A 


-0.479 706 380 368 904 997 8(2) 


-0.488 065 357 851 706 80(5) 


-0.496 762 894 249 561 327(1) 




-0.479 706 380 368 904 25 [2] 


-0.488 065 357 851 705 [2] 


-0.496 762 894 249 561 31 [2] 




-0.479 706 380 368 [IT] 


-0.488 065 357 841 [17] 


-0.496 762 894 248 \17} 



(a)'p]^g i^ggf; variational energies known from earlier calculations. 
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TABLE III: The total energies (E) of the bound S{L 
molecular ions in muon- atomic units (m^ = l,fi = l,e 
functions used in Eq.(2). 



= 0)— states in the symmetric muonic 
= 1). N designates the number of basis 



N 




E{dd^i) 


E{tt^£) 


3300 


-0.494386 820248 934808 76255 


-0.531111 135402 386448 65 


-0.546374 225613 816728 844 


3500 


-0.494386 820248 934808 76280 


-0.531111 135402 386449 61 


-0.546374 225613 816728 849 


3700 


-0.494386 820248 934808 76289 


-0.531111 135402 386450 59 


-0.546374 225613 816728 855 


3840 


-0.494386 820248 934808 76294 


-0.531111 135402 386451 22 


-0.546374 225613 816728 856 


A 


-0.494386 820248 934808 76325(10) 


-0.531111 135402 386455(2) 


-0.546374 225613 816728 90(3) 




-0.494386 820 248 934 546 94 [2] 


-0.531111 135402 386374 5 [2] 


-0.546 374 225 613 816 71 [2] 


E" 


-0.494386 820 248 58 \17} 


-0.531 111 135 402 [17] 


-0.546 374 225 598 [17] 



("^The best variational energies known from earlier calculations. 



TABLE IV: The total energies {E) of the bound P{L = 1)— states in the symmetric muonic 
molecular ions in muon-atomic units {m^ = 1, ?i = 1, e = 1). N designates the number of basis 
functions used in Eq.(2). 





E{ppfi) 


E{ddfi) 


Eittfi) 


3300 


-0.468 458 436 303 385 268 11 


-0.513 623 956 792 681 889 


-0.533 263 449 820 383 210 


3500 


-0.468 458 436 303 385 272 62 


-0.513 623 956 792 681 896 


-0.533 263 449 820 383 224 


3700 


-0.468 458 436 303 385 276 22 


-0.513 623 956 792 681 901 


-0.533 263 449 820 383 236 


3840 


-0.468 458 436 303 385 278 58 


-0.513 623 956 792 681 905 


-0.533 263 449 820 383 242 


A 


-0.468 458 436 303 385 35(5) 


-0.513 623 956 792 681 95(2) 


-0.533 263 449 820 383 28(2) 




-0.468 584 363 033 834 4 [2] 


-0.513 623 956 792 680 25 [2] 


-0.533 263 449 820 376 6 [2] 




-0.468 584 363 033 [17] 


-0.513 623 956 792 6 [17] 


-0.533 263 449 817 6 [17] 



('^^The best variational energies known from earlier calculations. 
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TABLE V: The total energies (E) of the bound D{L = 2)— states of the the dd^,ttfi and dtfi 
muonic molecular ions in muon atomic units. N designates the number of basis functions used in 
Eq.(2). 



N 


E{ddfi) 


E{ttii) 


N 


E{dtfi) 


3000 


-0.488 708 332 3489 


-0.512 568 653 4551 


3300 


-0.500 118 083 958 73 


3200 


-0.488 708 332 3559 


-0.512 568 653 5939 


3500 


-0.500 118 083 958 91 


3400 


-0.488 708 332 3604 


-0.512 568 653 7061 


3700 


-0.500 118 083 959 06 


3600 


-0.488 708 332 3625 


-0.512 568 653 7938 


3840 


-0.500 118 083 959 15 


3800 


-0.488 708 332 3638 


-0.512 568 653 8630 


4000 


-0.500 118 083 959 21 


A 


-0.488 708 332 375(10) 


-0.512 568 654 10(11) 


A 


-0.500 118 083 959 35(5) 




-0.488 708 332 3645 [2] 


-0.512 568 653 9435 [2j 


-0.500 118 083 959 201 [2] 




-0.488 708 331 28 [17] 


-0.500 118 083 90 [17] 


E" 


-0.512 568 651 99 [17] 



("^The best variational energies known from earlier calculations. 
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TABLE VI: The total variational energies in muon atomic units for all 22 bound states in the six 
muonic molecular ions. 



ion 


S{L = 0) -state 


S*{L = 0) -state 


ppfi 


-0.494 386 820 248 934 808 763 25(10) 




pdji 


-0.512 711 796 501 514 55(1) 




pt/i 


-0.519 880 089 782 921 0(5) 




dd/i 


-0.531 111 135 402 386 455(2) 


-0.479 706 380 368 904 997 8(2) 


dtfjL 


-0.538 594 975 061 413 3(2) 


-0.488 065 357 851 706 80(5) 


ttji 


-0.546 374 225 613 816 728 90(3) 


-0.496 762 894 249 561 327(1) 


ion 


P{L = 1) -state 


P*{L = l)-state 


ppfi 


-0.468 458 436 303 385 35(5) 




pdji 


-0.490 664 169 479 327(1) 




ptn 


-0.499 492 029 991 539(1) 




ddfi 


-0.513 623 956 792 681 95(2) 


-0.473 686 733 842 727 0(5) 


dtfj, 


-0.523 191 456 315 960(1) 


-0.481 991 529 973 85(5) 


ttji 


-0.533 263 449 820 383 28(2) 


-0.489 908 667 504 943 30(8) 


ion 


D{L = 2) -state 


F(L = 3) -state 


ddji 


-0.488 708 332 375(10) 




ddfjL 


-0.500 118 083 959 35(5) 




ttji 


-0.512 568 654 10(11) 


-0.490 554 165 50(25) 



TABLE VIL The convergence of the (r2i), (r'^i), {^21) and v^i expectation values for the ground 
(bound) S{L = 0)— states of some muonic molecular ions (in muon-atomic units). 





dtji 


pdfx 


ddfx 


N (r2i) 


(rl) 


{S32) 


1^31 


3300 3.299486184357381476285 


5.8818538929815795 


0.173456203086 


-0.9466714285 


3500 3.299486184357381476271 


5.8818538929815783 


0.173456202965 


-0.9466714315 


3700 3.299486184357381476272 


5.8818538929815769 


0.173456202754 


-0.9466714323 


3840 3.299486184357381476267 


5.8818538929815760 


0.173456202768 


-0.9466714337 



A 3.299486184357381476260(20) 5.8818538929815755(15) 0.173456202768(15) -0.9466714337(14) 



25 



TABLE VIII: The bound state properties X computed for the ground S{L = 0)— state and excited 
S*{L = 0)— state in the pd^ and dtji muonic molecular ions (in muon-atomic units). 



X 


pd^{S{L = 0)— state 


dtii(S{L = 0) -state 


dt4J.{S*{L = 0) -state 




0.6411463715926726(1) 


0.7227000085308685(3) 


0.5146887540974253(3) 


(^3-2') 


0.7533736132562925(1) 


0.7583156106645547(3) 


0.7053752991037914(3) 




0.3690963918459403(3) 


0.4038256690725974(5) 


0.2439333374978041(6) 


(rai) 


2.451487588757093(2) 


2.117912246542931(3) 


3.933235704904257(3) 




2.087699148755371(2) 


2.023720495653520(3) 


2.738751041467644(3) 


{r2l) 


3.100710403484118(3) 


2.747914131666301(5) 


5.161228965846044(5) 


(rl) 


8.033494173462087(3) 


5.881853892981577(4) 


22.39719316098318(5) 


(rl) 


11.082902111235402(4) 


8.2873253005241016(6) 


30.631300336084518(6) 


(rh) 


20.6547093169798(2) 


17.4696964524501(3) 


65.2551239488477(3) 


K) 


42.147963983631(4) 


27.208343544797(5) 


201.45178140015(6) 




0.2806191952124375(7) 


0.3910764762722844(9) 


0.3836327499261653(9) 




0.3674609039968256(7) 


0.4218837824381256(9) 


0.4476331914219544(9) 


0.46041133628342744(4) 0.50069529865293120(5) 0.50035761171633245(6) 


(5(r3i)> 


0.11770973301(5) 


0.15452554329(7) 


0.10723299252(7) 




-0.898787990(3) 


-0.946671420(5) 


-0.946671472(5) 


^31 


-0.8987879287819516 


-0.9466714310522288 


-0.9466714310522288 


(5(r32)> 


0.17345620275(5) 


0.17451466503(7) 


0.17875743059(7) 


^'32 


-0.9466716379(8) 


-0.9637483436(10) 


-0.9637483713(10) 


^32 


-0.9466714310522287 


-0.9637483334950246 


-0.9637483334950246 


(5(r2i)> 


1.4616938(4)-10-5 


8.870847(1)-10-^ 


7.416512(1)-10-^ 




5.9192044(3) 


10.644434(5) 


10.644863(5) 




5.9191833130515522 


10.644186704849556 


10.644186704849556 




2.2812()(2)-1() 


l.()l!)21(:-!)-l() *• 


l.;-!o81(i(:-!)-l() *' 



(")The exact value. 



26 



TABLE IX: Calculations of the mass gradients for the P{L = 1)— states in the ppfi,ddiJ, and tt^ 
muonic molecular ions in muon-atomic units where H = l,e = 1 and = 206.768262me = 1. The 
notations used in this Table correspond to the notations used in Eq.(6). 









ddfj, 




-2h) 


-0.46845 38220 26274 128170 


-0.51361 95148 47688 534232 


E{m^ 


-h) 


-0.46845 61291 64644 884357 


-0.51362 17147 27907 740353 


E{m^ 


+ h) 


-0.46846 07434 42495 278996 


-0.51362 61777 62223 263542 


E{m^ 


+ 2h) 


-0.46846 30505 81974 916597 


-0.51362 83987 29794 443654 


a 




-2.30713 89251 97(10)-10-3 


-2.23503 27015 23(10)-10-3 






//// 


(///')* 


EirUf, 


-2h) 


-0.53325 87089 79317 1418 


-0.48990 43312 78558 0219 


E{m^ 


-h) 


-0.53326 10794 00655 8570 


-0.48990 64993 92236 3403 


E{m^ 


+ h) 


-0.53326 58202 38499 2057 


-0.48991 08356 16741 7247 


E{m^ 


+ 2h) 


-0.53326 81906 55003 8459 


-0.48991 30037 27439 5589 


a 




-2.37041 89216 74(11)-10-^ 


-2.16811 22634 62(ll)-10-4 
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TABLE X: The total and binding energies of a weakly bound P*{L = 1)— state in the ddfi muonic 
molecular ion as the function of the muon mass rh^ = + y, where = 206.768262me and 
y is a non-negative real number. All energies have been recalculated to muon-atomic units where 
h = l,e = 1 and = 206.768262me = 1. The binding energies e are given in electrono Volts (eV). 
N is the total number of basis functions used in variational expansion. 



y 


E{N = 3300) 


e{N = 3300) 





-0.4736867338427258104 


-1.97498808799965 


1 


-0.4712335469779210548 


-1.68187156159084 


2 


-0.4688058829862271553 


-1.39988399201543 


3 


-0.4664034583422745329 


-1.12949080447425 


4 


-0.4640260218725323738 


-0.87131211085906 


5 


-0.4616733706509872075 


-0.62621871045415 


6 


-0.4593453849322646260 


-0.39554090739252 


7 


-0.4570421278654125672 


-0.18165985634448 


7.9 


-0.4549909488076416492 


-0.00875602540622 


7.95 


-0.4548776498295269995 


-0.00002551843452 


8 


-0.4547644307696878137 +0.00854711729854 


y 


E{N = 3840) 


e{N = 3840) 


7.95 


-0.4548776454620138526 


-0.00007865184962 
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TABLE XI: The total and binding energies of a weakly bound P*{L = 1)— state in the dtfi muonic 
molecular ion as the function of the muon mass rh^ = + y, where = 206.768262me and 
y is a non-negative real number. All energies have been recalculated to muon-atomic units where 
h = l,e = 1 and = 206.768262me = 1. The binding energies e are given in electrono Volts (eV). 
N is the total number of basis functions used in variational expansion. 



y E{N = 3700) e{N = 3700) 

0.0 -0.48199152997371597 -0.66033868534010 

1.0 -0.48184465428112405 -0.31071047997243 

1.8 -0.48173196328443396 -0.05573790640587 

1.9 -0.48171831993605128 -0.02623616728477 
1.95 -0.48171161233515334 -0.01212193523051 
1.98 -0.48170758889627615 -0.00365631831137 
2.0 -0.48170482239836190 +0.00246722246561 
y E{N = 3840) £{N = 3840) 

1.95 -0.48171161309408528 -0.01212624559609 
1.98 -0.48170759005686002 -0.00366291081502 
2.0 -0.48170492020870050 +0.00191157404752 
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